Optical amplification enhancement in photonic crystals 
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Improving and controlling the efficiency of a gain medium is one of the most challenging problems 
of laser research. By measuring the gain length in an opal based photonic crystal doped with laser 
dye, we demonstrate that optical amplification is more than twenty-fold enhanced along the F-K 
symmetry directions of the face centered cubic photonic crystal. These results are theoretically 
explained by directional variations of the density of states, providing a quantitative connection 
between density of the states and light amplification. 
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The study of unconventional lasing is a novel and 
active field of research. Mirror-less, micron-sized and 
three-dimensional lasers, in the form of random lasers [![, 
micro-sphere and microdrop lasers 0, photonic crystal 
(PhC) lasers Q and chaotic-cavity lasers [1] have been 
recently proposed as more sophisticated alternatives to 
standard lasers made of macroscopic mirrors and Fabry- 
Perot cavities. 

Artificially engineered photonic materials, with 
nanoscopic features and topologies ranging from order to 
disorder, allow for a superior control of light modes, dis- 
persion and gain [5j . Photonic structures can control and 
tune the directionality and spectral extent of the emit- 
ted light without relying on light absorption, allowing for 
efficient light sources with smarter functionalities. 

While the research of 2D PhC lasing can exploit high- 
quality and low volume micro-cavities @, in three di- 
mensions (3D) PhC lasing is so far only based on the 
large light-matter interaction that can be achieved by ex- 
tended Bloch modes, at band-edges 0] or at high-energy 
fiat bands @. These modes effectively confine the light 
in the active medium and increase the amplification and 
the gain extraction from large bulk volumes [9j. A simi- 
lar confinement, but originated by multiple scattering, is 
responsible for random lasing of diffusive modes [l] . The 
first pioneering study of gain in PhC has demonstrated 
a directional modification of the gain spectral profile, at- 
tributed to band-edge effects [l(| • Enhanced optical gain 
was also extrapolated from a decreased lasing threshold 
at the L-pseudogap ed ges of an opal-TiC^ composite con- 
taining quantum dots 

The lasing wavelength selection and the decrease of the 
lasing threshold in PhC lasers has been often claimed to 
be connected to spectral maxima of density of states [t| , 
but a quantitative evidence of the directional nature of 
the gain length due to the dependance of the density 
of the states on the crystallographic direction in which 
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FIG. 1: a) Photonic crystal transmission and reflection spec- 
tra. The shaded region indicates the gain curve of the dye. b) 
SEM image of the top of the sample with a sketch of the crys- 
tallographic directions and the real space arrangement. The 
hexagon indicates a cut of 3d Brillouin in the plane (111) c) 
SEM image of the side of the sample, corrected with the real 
color of the dye doped spheres. 
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the amplification grows has never been reported so far. 
Demonstrating this link is of paramount importance not 
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only for the future developments of nano-structured laser 
devices but also as a starting point for several fundamen- 
tal investigations ranging from quantum information to 
slow light phenomena. 

In this article we report a detailed measurement of 
the gain length in opal-based photonic crystals doped by 
Rhodamine 6G. The measured gain length has a six- fold 
symmetry, takes values up to ~ 500 cm^ 1 along the T- 
K directions in reciprocal space and decreases to values 
below measurable levels of ~ 20 cm -1 far from the high 
symmetry directions. We compare these findings to the 
directional density of states which has been predicted to 
be responsible for the increased gain. 

The samples under study are thin film artificial opals 
[lH consisting of an fee matrix of polystyrene spheres 
(PS, refractive index, n — 1.59, diameter d = 512 ± 15 
nm) bulk-doped with Rhodamine 6G with an emission 
centered around A = 574 nm (from Duke Scientific). Fig- 
ure [T] shows transmission and reflection spectrum from 
the sample (in the panel a) as well as a two SEM pic- 
tures, of the top of sample (in panel b) and of its side 
(in panel c). While a band-gap is visible at around 1200 
nm, the dye emission indicated by the gray area sits in 
the high energy peaks of the spectrum around 574 nm. 

The optical gain in the dye-doped three-dimensional 
artificial opals can be extracted by employing the stripe 
length technique [l3|, [l4{ . The emission intensity and am- 
plification is measured as a function of the length of the 
illuminated area on the sample. In figure [5^, a sketch 
of the measurement configurations used for the experi- 
ment is presented, while figure [^p report the curves from 
which the amount of optical amplification is extracted. 
The spontaneous emitted photons, traveling in the direc- 
tion z of the illuminated stripe, undergo an exponential 
amplification of the intensity with respect to the excita- 
tion length z as they pass trough an area of population 
inversion. This can be quantified by the gain coefficient, 
the most important actor in light amplification and las- 
ing, which is a robust quantity, independent on the sam- 
ple thickness (in absence of pump depletion) and on the 
edge coupling of the light from the stripe to the detector. 
The output intensity I{z) collected at one edge of the 
stripe is the sum of spontaneous and stimulated emission 
and has the well-known spatial dependence 



I(z) = I SE [exp(Gz)-l] 
G = g-l-^-i-\ 



(la) 
(lb) 



where Ise is the spontaneous emission (SE) collected by 
the detector, g is the gain (including confinement effects), 
£ a the absorption length and £ s the scattering mean free 
path. The quantity G, refereed to as the net gain, is di- 
rectly obtained from stripe-length measurements. Equa- 
tion [T] a, holds until gain saturation is reached, which 
occurs for stripe lengths larger than around 150-200 /im, 
in our case. 

In our experiment the excitation is carried out by illu- 
minating the sample along the PhC [111] direction (like 
in figure^) with a Q-switched Nd:Yag laser, 9 ns pulses, 




100 150 200 

Stripe lenght ( |jm) 

FIG. 2: (Color Online) a) Sketch of two configurations for 
measurement of the G parameter. The green laser pump is 
incident along the (111) direction and the red emission in 
the direction s(k, tJo) is collected on a side by a miniature 
spectrometer. The PhC is rotated continuously in the plane 
(111) orthogonal to the [111] direction, to explore the angle 
6 dependence of the emission. The illumination stripe has 
a variable length z. b) Intensity of the emitted light, as a 
function of the length, collected at different angles. Black 
squares correspond to the emission at 60 deg, blue triangles at 
246 deg, red triangles at 170 deg, while blue squares indicate 
the angle of 155 deg. The continuous lines represent the fits 
obtained by using the equation reported in the experimental 
section. 



10 Hz repetition rate, at a wavelength of 532 nm. The 
laser beam has the shape of a stripe of thickness ~ 20/im 
obtained by focusing the laser source with a cylindrical 
lens of 75 mm focal. Measurements of the dependence of 
J as a function of z has been performed by cutting the 
stripe with a blade moved automatically with micrometer 
accuracy. A miniature fiber-coupled spectrometer is used 
to efficiently collect the light along the stripe direction, 
as described in ref. [Til ]. 

By rotating the stripe relative to the in-plane direc- 
tion s, various crystallographic directions , including the 
high-symmetry ones, can be explored. In Fig. [5p, the 
dependence of the collected fluorescence intensity on the 
stripe length z, for various different orientations of the 
photonic crystal is shown. Emission for light propagat- 
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ing along two T-K crystal directions is presented (60 and 
240 deg), as well as an intermediate position (170 deg). 
A non symmetrical direction is probed too (155 deg). 
These measurements show a clear dependence of the flu- 
orescence intensity on the propagating angle in the crys- 
tal. 

This result is a consequence of the photonic crystal 
complex band structure that alter the optical gain and 
amplification processes as compared to an homogenous 
medium. In figure [3^, we report the isofrequency surface 
at the working frequency. For the sake of clarity only half 
the isofrequency, cut along a plane perpendicular to the 
[111] direction, is presented. Figure [3b shows the calcu- 
lated band structure along a path in the reciprocal space 
(sketched in panel a). From the full 3D band structure 
(calculated by preconditioned conjugate-gradient mini- 
mization of the block Rayleigh quotient in a planewave 
basis, using a freely available software package [HI) the 
isofrequency surface and group velocities are obtained. 




FIG. 3: (Color Online) a) Isofrequency surface section in 
the plane, for A = 574 nm, cut in the plane _L to the (111) 
direction, || to the sample surface, which is the plane of study. 
The color code refers to the band number indicated by n. b) 
Band diagram. The shaded red region indicates the gain curve 
of the dye centered at A = 574 nm, while the green line is the 
pump laser at A = 532 nm. Therelevant symmetry directions 
are shown by the capital letters. 

In the band diagram in figure [3Jd, the directions T— K, 
K-L and L'-r are explored when rotating the sample 
around the direction [111] normal to its surface. The 
shaded red region indicates the gain curve of the dye 
which excites the high order diffractive modes when 
pumped by a frequency doubled Nd:Yag laser (wave- 
length 532 nm) which is indicated by the green line. Note 



that U is not a high symmetry point equivalent to L, but 
just the midpoint of the experimental trajectory along 
the hexagonal facet of the Brillouin zone. 

An homogenous media would exhibit a circular (result- 
ing from the section of a sphere) isofrequency contour; in 
a photonic crystal a rich structure with a non-constant 
curvature is visible, which results from contributions of 
eight different bands around the emission energy. Each 
point in the isofrequency surface identify a wave-vector 
k, for which the direction of propagation of the energy s 
is given by the gradient of the surface. 
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FIG. 4: (Color Online) Measured gain in a PhC as a func- 
tion of the excitation angle. The orange thin line corresponds 
to the homogenous reference sample. The emission from the 
PhC is presented by the full symbols, while the gain of the 
reference sample is denoted by the empty orange circles. Dif- 
ferent symbol styles refer to measurements performed after a 
90 deg sample rotation. The curve above the gray shaded re- 
gion denotes the DOS. F-K directions occur each 60 deg from 
= deg. 

By fitting the data in figure [2Jd to equation [T] , the 
net gain G can be extracted and plotted as a function of 
the stripe direction as in figure 0] G is affected by the 
variations of group velocities and varies as a function of 
the probed crystallographic direction. The angle = 
corresponds to the configuration in which the stripe is 
aligned with the opal growth axis, and it identifies a T-K 
direction. A large variation of G is visible with a clear C6 
symmetry and maxima along the six T-K directions (each 
60 deg from 9 = deg) reaching values up to 500 cm -1 . 
The minima measured along less-symmetrical directions 
reach values which are limited to ~ 20 cm -1 by the mea- 
surement sensitivity: a more than 20-fold variation is ob- 
served. Different color symbols in figure 0] corresponds to 
measurement sets that have been performed by rotating 
the sample 90 degrees. Each set of data is thus obtained 
on a different side of the sample. The good quality of the 
latter - and, especially, the homogeneous dye distribution 
- allows a high reproducibility of the measurements, as 
shown in figure SI 

In order to prove that the observed effect depends on 
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the regular disposition of the spheres, a reference sample 
has been measured too: the same PhC sample after melt- 
ing by keeping at 100 Celsius overnight. This tempera- 
ture is enough to deform the PS spheres while keeping 
the film shape (as confirmed by SEM inspection) with- 
out degrading the dye molecules. The reference sample 
has a dye molecule density 35% higher than the regular 
sample, as it is a compact film as compared to the opal 
PhC whose filling fraction is 0.74. The measured gain of 
99 ± 20 cm -1 in the reference is now independent on the 
excitation angle. The orange thin line in figure [4] shows 
the reference sample gain multiplied by 0.74 to account 
for the dye molecule density increase. 

The data are compared to the directional density of 
states p(s,cj) that is calculated as the density of modes 
with Poynting vector pointing in the stripe direction, 
given by contribution of modes with any A:-vector in the 
3D space, and shown in figure 0] as the curve above the 
gray shaded region. 

The extraordinary increase in optical amplification and 
G in figure HI as predicted by Sakoda [§,[111], is connected 
to the variation of the density of states p (s, ui) in the di- 
rection s of the stripe. The number of available optical 
states in a given direction, i.e. for a given Poynting vec- 
tor parallel to s, determines the emission probability for 
the excited dye molecules. A higher value of p (s, uj) cor- 
responds to a higher emission probability in the (s,w) 
mode and therefore a larger gain for photons travelling 
in the s direction. The shaded region in figure 2] shows 
the variation of p (s, uj) along the various directions in the 
(111) plane. The theoretical curve in figure 0] has been 
rescaled by the ratio of the computed value of p (s, uj) 
for an homogenous medium to the measured gain in the 
reference melted opal. 

In general cases, also scattering and absorption can 
vary with the angle. Far from the dye absorption curve, 
the opal absorption length has been measured to be 
£ a ~ 10 m [17], which is therefore negligible for the sys- 
tem under study. The scattering mean free path £ s can 



present indeed large variations with p (s, uj), as studied in 
ref. [l~8l ]. For lower energy and similar opals it has been 
measured to be in the range 100 — 500 /im. This value 
is larger but comparable to the gain length £ g = G _1 
which is here in the range 20-200 ^m. Large scattering is 
anyhow responsible only for a decrease of the amplifica- 
tion and cannot explain the large increase in G that we 
have measured. If the variations of £ g were governed by 
i s , then the value of £ s should be smaller than £ g and no 
net gain could be measured. In our case scattering losses 
reduce the measured gain in figure [4] and set a minimum 
value of G that can be measured in our experiment (~20 
cm -1 ). 

In conclusion we have presented an observation of a 
six-fold increase of the gain in opal photonic crystal, as 
compared to the homogeneous film, and a more than 20- 
fold variation of G between T-K and less symmetrical 
directions, in the same photonic crystal. We explain this 
enhancement as due to a large increase of the density of 
the available modes when exciting around the T-K di- 
rection. Large variations of the gain in PhC show the 
impact of the tailored density of states on light gener- 
ation and amplification and open the way to enhance 
other phenomena like nonlinear wave mixing and har- 
monic generation. Our results show how nanostructured 
media could be at the basis of the development of novel 
lasing sources with exceptional tunability, directionality 
and efficiency while being plastic photonics CMOS com- 
patible, and candidates for in-board interconnections for 
future generation computers. 
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